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In inversion-asymmetric semiconductors, spin-orbit coupling induces a k-dependent spin splitting 
of valence and conduction bands, which is a well-known cause for spin decoherence in bulk and 
heterostructures. Manipulating nonequilibrium spin coherence in device applications thus requires 
understanding how valence and conduction band spin splitting affects carrier spin dynamics. This 
paper studies the relevance of this decoherence mechanism for collective intersubband spin-density 
excitations (SDEs) in quantum wells. A density-functional formalism for the linear spin-density 
matrix response is presented that describes SDEs in the conduction band of quantum wells with 
subbands that may be non-parabolic and spin-split due to bulk or structural inversion asymmetry 
(Rashba effect). As an example, we consider a 40 nm GaAs/Alo.3Gao.7As quantum well, including 
Rashba spin splitting of the conduction subbands. We find a coupling and wavevector-dependent 
splitting of the longitudinal and transverse SDEs. However, decoherence of the SDEs is not de- 
termined by subband spin splitting, due to collective effects arising from dynamical exchange and 
correlation. 

PACS numbers: 71.15.Mb; 71.45.Gm; 73.21.Fg; 72.25.Rb 



I. INTRODUCTION 

Most currently available semiconductor device tech- 
nologies are entirely based on manipulating electronic 
charges. The emerging field of spintronicsfclcl proposes to 
exploit, in addition, the spin degree of freedom of carri- 
ers, thereby adding new features and functionalities to 
solid-state devices. Many of the proposed new appli- 
cations rely, in one form or another, on manipulating 
nonequilibrium spin coherence. The hope that this may 
indeed lead to viable practical approaches has beeij_aiip- 
ported through recent experimental observationsEruCI of 
long-lived (> 100 ns) and spatially extended (> 100 /im) 
coherent spin states in semiconductors. Two character- 
istic times, T\ and T2, provide a quantitative measure 
for the magnitude and persistence of spin coherence. T\ 
describes the return to equilibrium of a non-equilibrium 
spin population, and T2 measures the coherence loss due 
to dephasing of transverse spin order (for more details, 
see Ref. |). 

Spin relaxation in GaAs quantum., . we lls was 
recently studied experimentallyul§BE3llll and 
theoretically.E3E3o'Ej'E£l Measurements of the elec- 
tronic T\ involve circularly polarized pump-probe 
techniques to create and observe coherent spin pop- 
ulations in the lowest conduction subband. Electron 
spin decoherence has been shown to occur via the spin 
precessionEH of carriers with finite crystal momentum k 
in the effective k-dependent crsiatal magnetic field of an 
inversion-asymmetric matcrial.Ea We note that the the- 
ory of Refs. |^,[l5| gives good agreement with experiment, 
without including electron-electron interactions. 

In this paper, we will consider electronic charge and 
spin dynamics in quantum wells involving not one but 



two subbands (we will limit the discussion here mainly 
to conduction subbands). One motivation for this work 
is that intersubband (ISB) charge dynamics in quantum 
wells is .currently of great experimental and theoretical 
interest Jl3 since electronic ISB transitions are the basis 
of a variety of new devices operatise in the terphertz 
frequency regime, such as-detectors,E2l modulators^ and 
quantum cascade lasers. E3o In view of this, it seems 
worthwhile to explore ISB spin dynamics as a possible 
route towards novel applications in the terahertz regime. 

Analogous to the case of spin dynamics discussed 
above, one may define characteristic times for intersub- 
band dynamics (for an overview, see Ref. p4| ). Population 
decay from an excited to a lower conduction subband 
is measured by an intersubband relaxation time T^ 5 , 
and loss of coherence of collective intersubband excita- 
tions is measured by a dephasing time T^ 13 . These two 
times have been measured experiment ally f or ISB charge- 
density excitations in quantum wellsc^E3 and found to 
differ substantially at low temperatures, being three 
orders of magnitude smaller than T-[ SB . The reason is 
that intersubband relaxation proceeds mainly via phonon 
emission and is thus slowed down by an energy bottle- 
neck for small phonon momenta as well as for the op- 
tical phonon branch. This differs from the case of con- 
duction electron spin relaxation, where T\ and T2 are 
comparable.EftLj 

On the other hand, dephasing of collective ISB exci- 
tations in quantum wells is determined by a complex in- 
terplay of a variety of different scattering mechanisms, 
whose relative importance is not a priori obvious. In re- 
cent experimental^ and theoreticalEa work, it was found 
that the linewidth of (homogeneously broadened) ISB 
charge plasmons in a wide GaAsZAlo.3Gao.7As quan- 
tum well, where phonon scattering plays no role, is de- 
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termined mainly by interface roughness and electronic 
many-body effects. 

The question now arises which physical mechanisms 
govern the dephasing of collective ISB spin-density exci- 
tations. As a first step towards a clarification of this 
question, this paper addresses the influence of the In- 
dependent crystal magnetic field in semiconductor quan- 
tum wells on ISB spin-density excitations, and the impor- 
tance of many-body effects. The latter will be described 
in the framework of density- functional theory (DFT). 

This paper is organized as follows. In Section II 
we set up the formalism for calculating the electronic 
ground state in modulation-doped quantum well conduc- 
tion subbands, includingspin-orbit coupling and many- 
body effects. Section III presents a general response 
formalism for the spin-density matrix, based on time- 



dependent DFT (TDDFT). In Section [TV] we consider an 
explicit example and calculate the collective ISB charge- 
and spin-density excitations in the conduction band of 
a GaAs/Alo.3Gao.7As quantum well, including spin-orbit 
coupling. Section ^ contains our conclusions. Various 
technical details can be found in Appendices [a| and ||. 



II. ELECTRONIC GROUND STATE 

We consider a modulation-doped quantum well (direc- 
tion of growth: z-axis) containing N conduction electrons 
per unit area. In th e s tandar d multi-band k ■ p approach 
for semiconductors, E^EffillH the single-particle states in 
a quantum well are expanded in terms of Bloch functions 
at the zone center, u n (r): 



*iq„( r ) 



tq,|i 



4>jn(z) u n (r) 



(1) 



where ipj n (z) are envelope function belonging to the jth 
subband, and rn = (x,y) and q|| = (q Xl q y ) are in-plane 
position and wave vectors. In general, u n (r) =u n ^(r)^+ 
u n i (r)£j , where denote two-component Pauli spinors. 
The functions u na (r) have the form u na (r) = c na \j) 1 
where c na are complex coefficients. Usually Nb includes 
at least the conduction-band *-states and the valence- 
band p states (8- band or KanecJ model), but in general 
a 14-band model is needed for a cpasistent description of 
spin splitting in heterostructures.tSl In that case, the set 
of basis functions is 

|i) 6 [\Z) v ,\X±iY) v ,\S) c ,\Z) c ,\X±iY) c ] . (2) 

This leads to a Hamiltonian in 8 x 8 (pti-ljl X 14) matrix 
form, whose elements are well known£3H The envelope 
functions ipj n (z) for valence and conduction bands then 
follow frorpjthc resulting 8 (or 14) coupled single-particle 
equations .EJ 

If one is only interested in the electronic structure of 
the conduction band of a quantum well, it is convenient 
to reduce the multi-band Hamiltonian described above 



to a 2 x 2 conduction band Hamiltonian.@0Si»O The 
single-particle states (|l|) can then be simplified to the 
following two-component form: 



(3) 



The envelope functions <pj a follow from a two-component 
effective-mass Kohn-Sham equation: 



E 

13=1,1 



h S af3 + v e ^(z) + H s ° p {z) + v^{z) ^(q,| , z) 



Ejq u <Pj a (q\\,z) 



(4) 



h 2 qf. 



where a =|, |, and 

d ft d 

h = -d~z tmeJ) Tz + MM + Vcont {z) + m{z) - 

(5) 

The spin-independent part h of the 2x2 conduction 
band Hamiltonian accounts for possible non-parabolicity 
of the subbands through an effective mass that depends 
on Ejq u . Explicit expressions for m(E, z) can be found in 



Refs. 36:37]. v con f(z) is the confining bare quantum well 
potential (e.g., a square well), and the Hartree potential 
Vn(z) is related to the electron ground-state density n(z), 
defined below, through Poisson's equation: 

d 2 v u {z) 



dz 2 



-Aire* n(z) , 



(6) 



where e* = e/^/e is the effective charge (e is the static 
dielectric constant of the material). 

Let us now discuss the spin-dependent parts of the 
Hamiltonian in Eq. (||). The first term, v^(z), describes 
externally applied uniform static electric and magnetic 
fields, E and B: 



eE z z S a p + i g*{z)fi B B • a , (7) 



where E = e z E z is perpendicular to the quantum well, 
and B can have arbitrary direction, a is the vector of 
the Pauli spin matrices, and g*(z) denotes the g- factor 
of the bulk material at point z. 

Intrinsic conduction band spin splitting, caused by 
spin-orbit interaction, in general comes from several dif- 
ferent sources. One often deals with situation where 
there are two major contributions: = + H^p, 
where BIA and SIA denote bulk and structural inver- 
sion asymmetry. The first term has the well-known form 

= [Mi • a/2] a0 , where ft = j(q x (ql - q 2 z ),q v (q 2 z - 



EJBIA 

H af3 



for bulk zincblende semiconductors 



E 



For a quantum well, H™^ depends on the growth direc- 
tion. For instance, along [001] we have 



H*\ A = [ H 



Ml 



1 d-y 
2~dz 



1 2 y) ( ^37+737 



d_ 

dz 



ETBIA 



= H 



Ml 



d d 

dz dz 



+ iq y ) 



(8) 
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for [110] and [111] directions, see Ref. [35| The second 
contribution to intrinsic spin splitting, SIA (also known 
as the Rashba effectES), has the form 



TTSIA 
rrSIA 



rrSlA 







i drj 
"2dz 



{Qx - iq v ) 



(9) 
(10) 



The material parameters 7(2) and T](z) are explicitly 
given in Ref. ^57]. We also mention a possible additional 
source of spin splitting in quantum wells, the so-called 
native interface asymmetry, related to chemical bonds 
across interfaces (for more details see Refs. |s|j39|). 

A novel feature in Eq. (Q), which distinguishes the 
present approach from previous studies of conduction 
band non-parabolicity and spin splitting, is that many- 
body effects are explicitly included through the exchange- 
correlation (xc) potential v* c g (z). xc effects have pre- 
viously been shown to produce non-negligible shifts of 
quantum well .subband energies and ISB charge plasmon 
frequencies.oESa As discussed below, including xc ef- 
fects is crucial for a physically correct description of col- 
lective ISB spin excitations. 

The solutions of Eq. (|j) have the interesting prop- 
erty of being mixed spin-up and spin-down eigenstates, 
due to the off-diagonal terms in the Hamiltonian caused 
by spin-orbit coupling and, possibly, externally applied 
transverse magnetic fields. The off-diagonal terms in H^g 
depend on q| | , and there is no choice of basis which di- 
agonalizes Eq. (|4|) for all q| | . Due to the absence of a 
global quantization axis, spin is no longer a good quan- 
tum number. This, requires a generalization of the well- 
known spin-DFTcj to systems with non-collinear spin. 
So far, this was done at only few occasions in the litera- 
ture, namely for non-collmeax magnetic materials such as 
7-Fe, U3Pt4 and Mft 3 Sn,E3a and inhomogeneous quan- 
tum Hall systems,C3 but, to our knowledge, never before 
in the present context of semiconductor nanostructures. 

Formally, the xc potential is defined as 



"a/3 



(r) 



SEL 



Sn a/ 3{r) 



(11) 



where the xc energy of the-system, E xc [n], is a functional 
of the spin-density matrixes 



a(r) = ££in*iiii*iq,i = 



Iff Tin 



(12) 



where fjq.. = f(Ep—Ejq..) denotes the Fermi occupation 
function, and Ep is the conduction band Fermi level. For 
*j q|| given by Eq. (0), we have 



Ml 



( Z ) = 2/iqi|l^t(<l||>*)| S 



[similarly for and 



n u( z ) 



(*) = ^/ J - q|l ^ T (qih2)^i(qn^) ■ (14) 



The usual approximation is 
( |l2| ) to be locally diagonal, 
non-collinear spin reads 



ake the density matrix 
44 so that the LSDA for 



On 



af3 



(15) 



3 hom 



(n, £) is the xc energy per particle of a homogeneous 
electron gas of density n and spin polarization £, which, 
is well known from quantum Monte Carlo calculations .Ej 
The local density and spin polarization are given by 



Tr n 



£ = — Tr an 
n 



(16) 
(17) 



The ground-state density is normalized as / dz n(z) — 
N . Explicit expressions for i)™ B (z) are given in Appendix 
[A|. With this form for v™ g {z), the 2x2 effective-mass 
Kohn-Sham equation (|J) is now completely defined. Self- 
consistent solution yields a set of subbands which are 
occupied up to Ep. 



III. LINEAR RESPONSE FORMALISM FOR 
THE SPIN-DENSITY MATRIX 

Once the electronic ground state (characterized by a 
set of subband levels and wave functions) has been calcu- 
lated, the next step is to consider excitations. The formal 
framework for describing excitations in electronic ra&iau 
body systems is provided by linear response theoryE3^H 

For the case where the wave functions take on a two- 
component form, the TDDFT linear response equation 
for quantum wells becomes a 2 x 2 matrix equation: 



™™<( k |h^ w ) = I dz> ' Xal',\\'( k \\> z > z ' ' w ) 

A,A' = T,i J 

xi$(k|,,*', W ). (18) 



This expresses, formally exactly, the first-order change of 
the spin-density matrix, n^, , via the response of a non- 
interacting system, characterized by the response func- 
tion Xaa'.w ( see below), to an effective perturbing po- 



tential of the form v 



(l.cxt) 

,J w> 



(1) 

XX' 



V 



(l,ext) , tt,H) 



XX' 



XX' 



(l,xc) 

'Ax' ■ 



Here, 



is the external perturbation, and the linearized 



(13) Hartree and xc potentials are 



4 



(l,Hl/, / \ (l,xc)/, / \ 
' J \\> ( k |h z > w ) + v \\> ( k |h z 



E 

C,C'=T,I 



dz" 



2ne* 



-k u \z'-z"\ 



(19) 



In the widely used adiabatic local-density approximation (ALDA).E3 the xc kernel is given by 

Q2 homf \A\ 

f&, cc ,(k N , z, z', u) = ° <K* - • 

The non-interacting response function takes on the form of a fourth-rank tensor: 

X Kb (r,r» = ^ E (/ill l - /iqf,) " 



(20) 



(21) 



where the *f?'s are given by Eq. ([!]). Formally, this multiband response formalism describes transitions among 
valence and conduction subbands, as well as interband transitions. In this paper, however, we focus exclusively on 
intersubband transitions in the conduction band of modulation doped heterostructures. Using Eq. (^), one can then 
transform the response function (Ell) into 



Xaa',AA'( k |b2>2'. W ) = E / 7^ 

jl J [ 



d 2< l\\ /iq||-k|| /jq|| 

(2tt) 2 u - E jcL . + E icl , ,_ k .. + ir] 



xip ja (q lll z)ip^,(q ll - k||,z)^* A (q||,z')^A'(q|| -k||,z') , 
I 



(22) 



where the Kohn-Sham envelope functions ipj a and ener- 
gies Ejq.. are obtained from Eq. (§). 

One can combine the perturbing spin-dependent po- 
tentials v„L and the solutions tllI, of the response equa- 



tion (18) in the fo llow ing, p hys ically more transparent 
way, see also Eqs. (A3) and (A4) of Appendix [A|: 



Vf> = Tr 



'3 iL 



(23) 



In terms of these quantities, the linear response equa- 
tion (|Tq) takes on the following form: 



3 

mf\k ll ,z,u) = J2 /dz'n^(k||,z,z',c)U fc (1) (k||,z',c). 

(25) 

The response functions 11^ and Xaa< w are related as 
follows: 



(1) rp 

uij = Ir 



(24) 



j = 0,1,2,3, where erg is the 2x2 unit matrix, and 
(j i , (72, (73 are the Pauli matrices. 



,(i),„(i) 



describes a collective charge-density 



excitation (CDE), and m 



(i) _ (i) 
- n n 



Ml 



is a longitu- 



dinal spin-density excitation (SDE) with respect to the 
z-axis. In terms of this choice of global spin quantiza- 
tion, mfp = ni 1 ) + njl and = i[n 



(i) 
11 



'It- 



ap- 



pear as transverse spin-density (or spin-flip) excitations. 
The CDE couples to an oscillating electric field polarized 
along the ^-direction, associated with Vq . The longi- 
tudinal SDE is excited by an oscillating magnetic field 
along z associated with V3 , and the transverse SDEs 
are excited by oscillating magnetic fields along x and y, 
associated with and V% , respectively. We will dis- 
cuss these selection rules in more detail below. 



TTKS _ KS 



dm 



(i)' 



(26) 



arc 



where the coefficients dm^ / dn a l, = dVj /dv^J., 
easily obtained from Eqs. ( p3| ) or (23). Explicit expres- 
sions for are given in Appendix B| The V% , in turn, 
are given as sums of external perturbations and linearized 
Hartree and xc terms: 



V fe (1) (k||,z, W ) = V fe cxt (k|h^) 



1=0 



EK 



2ire* 



kO°lO 



/«(k||,2,2',w) 



m, (1) (k||,z',w). (27) 



The xc kernels in ALDA are given in Appendix 



■5 



IV. RESULTS AND DISCUSSION 
A. Kohn-Sham wavefunctions and Rashba effect 

We will now discuss an example to illustrate the 
spin-density matrix response formalism developed above. 
Consider the case ofr-a,r4£) nm wide GaAs/Alo.3Gao.7As 
square quantum well,r 8 tr°l without any externally applied 
static electric and magnetic fields. We make the simpli- 
fying assumption of parabolic conduction subbands (i.e., 
neglecting the difference of the effective masses in well 
and barriers). Furthermore, we neglect B1A, but assume 
spin splitting is dominated by SIA, described by a sim- 
plified Rashba term of the formEa 



H s 



a[a x qj z 



R 

R* 



(28) 



where R = a(q y + iq x ), and a is taken to be a real, pos- 
itive adjustable parameter. The Rashba field is thus as- 
sumed to be the same for all conduction subbands, which 
is a reasonable approximation for wide quantum wells. 
The two-component Kohn-Sham equation (H) becomes 




E. 



Ipil 



where i = 1, 2, 3, . . ., and 

1 / d 2 

dz 2 



h 



2m* 



—-2 + q u + Uconf (z) + vh(z) 



Equation 



is solved by the following ansatz: 



Vvrfaib 2 ) 



1 



V2 



<Pj(z) 



s R* 
V5W\ 



tp.j{z) , 



(29) 



(30) 



(31) 



(32) 



where we replaced the subband index i by the pair of 
indices {sj}, such that s = (— l) 1 and j = (i + l)/2 for 
i odd and j — i/2 for i even. In the absence of the off- 
diagonal terms in Eq. (p9|), i.e. for inversion symmetry 
and hence spin degeneracy at each q||, j simply labels 
the spin-degenerate pairs, and s labels the members of 
the pairs. 

It is not difficult to see that in the presence of R the 
ground-state density matrix remains diagonal with = 



nn, an d hence = uf? = ^xc and v?f — trf? = 0. 



"TT — II — xc U ~ IT 

The fj{z) are therefore simply the solutions of the spin- 

unpolarized effective-mass Kohn-Sham equation 



1 d 2 

2m* dz 2 



(33) 



where ej are the energy levels of the associated, doubly 
degenerate, parabolic subbands. The presence of the off- 
diagonal Rashba terms in Eq. (|29|), however, lifts the spin 
degeneracy for q|| 7^ 0. We thus obtain, using \R\ = aqu, 



E. 



2m* 



saq\\ 



s = ±1 



(34) 
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FIG. 1: Lowest (1— »2) ISB excitation frequencies versus elec- 
tronic sheet density A in a 40 nm GaAs/Alo.3Gao.7As quan- 
tum well, at k|| = 0, for values of N where only the low- 
est subband is occupied (with parabolic subbands and with- 
out spin-orbit splitting). Full line: single-particle excitations 
(uj — €2 — ei). Dotted line: charge-density excitations in RPA. 
Dash-dotted and dashed lines: charge- and spin-density exci- 
tations in ALDA. 



for the energy eigenvalues associated with the solutions 
(HD^ofEq. |D. 



B. Collective intersubband excitations 

In the following, we will consider only cases where the 
lowest conduction band is occupied, which restricts the 
electron density in the quantum well to A < 1.82 x 
10 11 cm -2 . The goal is to study collective charge- and 
spin-density excitations between the first and the second 
subband. These collective modes are obtained by solv- 
ing the response equation (25) for the case where the 
external perturbation is zero. In that case, the effective 



perturbing potential consists of the self-consistent 
linearized Hartree and xc terms only. 

Let us first consider the case without spin-orbit cou- 
pling (a = 0). Figure |l| shows the density dependence 
of various ISB excitations at zero in-plane wavevector 
(k|| = 0). The full line depicts the single-particle excita- 
tions with frequencies w = £2 — ei, i.e. the bare Kohn- 
Sham excitation energies. The dotted line shows the ISB 
charge-density excitation in RPA, i.e. setting fgf = in 



the effective potential Eq. (|27|). The RPA excita- 

tion energies are always higher than the single-particle 
excitations, due to the so-called depolarization shift.EJ 
The ISB charge-density excitation in ALDA is shown by 
the dash-dotted line. Including xc effects in the response 
calculation produces a downshift of the plasmon energy 
of up to 0.75 meV. Finally, the spin-density ISB excita- 
tion is shown by the dashed line. In RPA, this excitation 
coincides with the single-particle excitation, since the de- 
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FIG. 2: ISB charge and spin plasmon wavevector dispersions 
in a 40 nm GaAs/Alo.3Gao.7As quantum well, for Rashba co- 
efficients a = (top) and a = 10 meVA (bottom). The 
electronic sheet density is I x I0 11 cm -2 . The shaded regions 
indicate Landau damping of the charge and spin plasmons. 
For a — 0, both regions coincide. For finite a, the Landau 
damping region for charge plasmons is unchanged (darker re- 
gion), but grows for spin plasmons (darker plus lighter region). 
The charge plasmon is essentially independent of a, but the 
spin plasmon splits into three branches for finite a. 



polarization shift affects only the charge mode. Thus, the 
spin plasmon only exists as a distinct, collective excita- 
tion because of xc effects. 

We now include spin-orbit coupling in the quantum 
well material by taking a finite, density-independent 
value of a = 10 meVA for the Rashba coupling param- 
eter. This is a typical value for practical situations of 
interest, for instance when applying a static electric field 
of strength 10 kV/cm in a GaAs quantum wellEa 

The ALDA in-plane wavevector dispersions of the ISB 
plasmons are shown in Figure |^, comparing the case of 
a = (top) and finite a (bottom). The shaded regions 
indicate Landau damping, i.e., collective modes overlap 
with the particle-hole continuum and can decay into in- 
coherent particle-hole pairs. In both cases, the charge 
plasmon lies above the region of Landau damping, and 
the spin plasmons lie below. In the case of a = 0, 
there is a common region of Landau damping for the 




-0.15 



0.0005 0.001 

fc„ (A)- 1 



0.0015 



FIG. 3: Splitting of the ISB spinplasmon dispersions, for the 
same quantum well as in Figure 0. Alu denotes the difference 
of the spin plasmon frequencies at a — 10 meVA and a — 
0. The dots indicate that the plasmons enter the region of 
Landau damping. The inset illustrates the selection rules (see 



text). 



.(1) 



and wjy 



are coupled and two-fold split. To lowest 



order in a, the splitting has the form S = C(N)aku, where 
C(N) depends on the electron sheet density. 



charge- and the spin plasmons. For finite a, the region 
of Landau damping for the spin plasmons grows, while 
for the charge plasmons it stays unchanged. In the ab- 
sence of other intrinsic or extrinsic scattering mechanisms 
(phonons, disorder), all collective modes outside the re- 
gion of Landau damping have infinite lifetime in ALDA. 

The charge plasmon dispersion is essentially indepen- 
dent of a. The spin plasmon, however, splits up into 
three branches for finite a. This is shown in more detail 
in Figure ^, where Alv denotes the difference of the spin 
plasmon frequencies at a = 10 meVA and a = 0. There 
are three different spin plasmon modes, all degenerate at 
a = 0. We will now discuss the nature of these modes, 
and how they couple to external fields. 

The charge and spin plasmons with k|| =0 cou- 
ple to external spin-dependent potentials of the form 
<A cxt )(z,w) = eEoza-j. For j = Q (CDE), y(Lext) is re _ 
lated to an oscillating uniform electric field Eexp(— iuit), 
where E = E e z (linearly polarized along z, perpen- 
dicular to the quantum well plane). For the SDE's 
(j = 1,2,3), ?/ 1,cxt ) corresponds by comparison with 
Eq. (Q) to oscillating magnetic fields B exp(— iuit), where 
B = 2E z/g*(z)nB- The CDE and SDE's can thus 
be formally viewed as collective electric and magnetic 
dipole transitions. At finite kn. the plasmons couple 
to external potentials of the formEJ uC 1 ' 6 **) (kn , z, to) = 
ei?o a o exp(kuz)aj, where <Zg is the effective Bohr radius. 
As before, these potentials can be related to oscillating 
electric and magnetic fields, Eexp[z(k|| • r|| — Lot)] and 
B exp[i(k|j • r|| — u)t)\. 

The inset in Figure || illustrates the selection rules for 
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TV (10 11 cm" 2 ) 

FIG. 4: ISB spin plasmon splitting coefficient C(N) (see 
Figure ^ , versus electron sheet density. 



the individual SDE modes: (i) a longitudinal mode, de- 
noted by mi , which couples to a magnetic field per- 
pendicular to the quantum well, B = B z e z . (ii) two 

transverse (or spin-flip) modes, mfj^ and , which 
couple to magnetic fields in the plane of the quantum 
well, B = Bi^ee and B = B\\ t £t, where £i = kii/fcii and 

£ t = £( x £ z . Figure [5] shows that, at finite a, mi and 



in 



(i) 
\\f 
(i) 
'II' 



are coupled and two-fold split. On the other hand, 

depends only very little on a, except a small red- 
shift independent of fcn. This small redshift, as well as 
the small splitting between the z and | \t modes at fc|| — 0, 
can be shown to be proportional to a 2 . 

Writing the Hamiltonian ( p8| ) in the form H S1A = 
hflji ■ a/2 defines the Rashba effective magnetic field 
£Ir = (2a/h)(q y , — (? x ,0), which lies in the quantum well 
plane and is perpendicular to qi i . Since in our example 
all subbands experience the same Hr, a collective ISB 
excitation with wavevector k|| implies a change in the 
effective magnetic field A£Ir = (2a/h)(k y ,—k x ,0) for 
all single-particle transitions, where Ai~2/e | |e* (see Figure 
||). This explains the physical origin of the splitting be- 
tween the different SDE branches: The two spin plasmon 
branches whose energies are shifted {z and | \t) are those 
responding to fields perpendicular to AH^, whereas the 
one which to lowest order in a does not shift (\\t) is par- 
allel to AJ7jj. Thus a spin polarization in either z or 
| \£ will precess in the z — £ plane. There are two possi- 
ble linear combinations, mi 1 ' ± im^J, one precessing in 
that direction which is favored by AJljj, the other in the 
opposite direction, thus costing more energy. 

To lowest order in a, the magnitude of the splitting 
between the two linear combinations of mi and mf,^ , S, 
is proportional to a and grows linearly with fen, i.e. S = 
C(N)ak\\. The numerical prefactor C(N) is a function 
of the electron density, and is plotted in Figure ||. The 
SDE splitting strongly increases with electron density, 



and reaches values of around 0.1 — 0.2 meV for sheet 
densities of order 1 x 10 11 cm -2 and higher, and plasmon 
wavevectors of order 0.001 A~ x . These splittings should 
be experimentally observable ..for example using inelastic 
light scattering techniques l 51 ^ 2 ! This would provide an 
opportunity for measuring the Rashba coefficient a. 



CONCLUSION 



We have presented a microscopic theory of collec- 
tive charge- and spin-density excitations in semiconduc- 
tor quantum wells based on spin-density-functional the- 
ory, with specific emphasis on intersubband excitations 
within the conduction band. The approach consists 
of two steps. We first calculate the electronic ground 
state in the quantum well (subband levels and envelope 
functions), including conduction band non-parabolicity 
and spin splitting, which leads to a 2 x 2 conduction 
band Hamiltonian. The associated Kohn-Sham matrix 
equation features spin-dependent xc potentials which are 
functionals of the spin-density matrix. 

We then determine the excitation energies using linear 
response theory in the formulation of TDDFT. Formally, 
one needs to solve a 2 x 2 matrix equation for the cou- 
pled charge- and spin-density-matrix response, including 
dynamic many-body effects. 

To illustrate the formalism, we considered the exam- 
ple of a quantum well with parabolic subbands that are 
split by a Rashba effective magnetic field. The charge 
plasmons were found to be independent of the Rashba 
field. The three possible spin plasmons, which are degen- 
erate in the absence of spin-orbit coupling, were found to 
be split into three branches, the splitting being propor- 
tional to the in-plane wavevector and to the strength of 
the Rashba field. 

This study illustrates the importance of including 
many-body effects beyond the RPA in calculating col- 
lective spin excitations. The collective nature of the 
ISB spin plasmons is purely a consequence of dynam- 
ical xc effects. Due to these collective effects, is 
not influenced by the precessional dccohcrcncc mecha- 
nisms related to spin-orbit coupling-ssjhich determine the 
intraband spin relaxation time T2.EZI Therefore, in the 
absence of impurities, disorder and phonon scattering, 
the lifetime of ISB spin plasmons is limited by dynami- 
cal many-body effects only. To capture these effects, one 
has to go beyond the ALDA and include retardation.Ej 
Within the ALDA, on the other hand, collective CDE 
and SDE's are infinitely long-lived. 

The effect of non-parabolic bands and more general 
forms of spin-orbit splitting (both BIA and SIA) in semi- 
conductor quantum wells will be addressed in more quan- 
titative detail in the future. 
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</3 = " 



finds 



i=0 



dn a 



1 „,xc 
i ) 



a,(3 =1,1 , 



(A9) 



APPENDIX A: LSDA FOR NONCOLLINEAR 
SPINS 

In LSDA0 the xc energy per particle of a spin- 
polarized homogeneous electron gas, e^ c is usu- 
ally appjpximated by the von Barth-Hedin parametri- 
zationBE-l 



<9e x 



TT 



7 U 



dn 



de x 



On 



l TT ~ n ll 



+ c 



de xc 
dC 



d( 



(A10) 



(All) 



eZ c (n,O=et(n,0) + [et(n,l)-et(n,0)}m)- (Al) 



The interpolation function between the paramagnetic 
(£ = 0) and ferromagnetic (£ = 1) limits, 



/(0 



(l + 4 / 3 + (l-0 4/3 

2(2 1 /3_i) 



(A2) 



reproduces the exact ^-dependence of the exchange en- 
ergy of the homogeneous electron gas, and approximates 
the ^-dependence of the correlation enetgj!,.— ..— . 

In the LSDA for non-collinear spins Jld'c3'L§C3 one still 
uses the spin-polarized homogeneous electron gas as ref- 
erence system, assuming that the xc energy per particle 
depends only on the local ground-state density n and the 
absolute value of the local ground-state spin polarization 
£, where, using definitions ([r| and (|l|), 



n = n TT + nu 



(A3) 



J n 



i r 



2n iT de 



n( d( 



2n n de xc 



n( d( 



(A12) 



(A13) 



Equations ( AlO)-.dAl3) are in agreement with the results 
of Heinonen et o,Ie3 

Next, we calculate, in ALDA, the xc kernels needed in 
Equation (p7|). The definition is 



d 2 [ne* c (n,Q] 
drrijdmk 



mi-mi(z) 
i=0,l,2,3 



S(z - z') , 

(A14) 



with the following results [omitting the 8{z — z')]: 



rxc 
100 



de h 
dn 



dn 2 dndt; n 9£ 2 



(A15) 



«T1 + n lT 

i{nn - "it) 




so that 



rn 



3 i 



(A4) 



(A5) 



and the sign of £ is determined with respect to the chosen 
global quantization axis: sign(£) = sign(rif| — n^). 

The xc potential in LSDA may then be obtained as 
follows: Letting n = mo, we define 



drrij 



mi —mi(z) 
1=0,1,2,3 



which yields 



det 



det 



i = 1,2,3. 



(A6) 

(A7) 
(A8) 



rxc 

JOi — 



mi d 2 e>> c 
n£ dndt; 



rXC 

J ij 



8^ de* 



niinij 
< ~dT ~ «) 3 



(He 



de 



(A16) 



(A17) 



Jji 



where i,j = 1,2,3 in ( |A16|) and (|A17p , and /* c 
for all i,j. For spin unpolarized ground states (mi 
m,2 = = 0), only those xc kernels diagonal in i, j are 
nonzero, with 



rxc _ r, 

J oo — - 



3e£ c (n,0) 



On 



d 2 e* c (n,0) 
dn 2 



and 



rXC 

J a 



4/9 



i(2V3 - l) 



[e£ c (n,l)-e£ c (n,0)] 



(A18) 



(A19) 



for i = 1, 2, 3. Notice that /qq and fff have the same ex- 
change parts, / X = fl = (4/9n)e£(n,0) = -(Qirn 2 )- 1 / 3 , 
but in general have different correlation parts. 
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APPENDIX B: NON-INTERACTING RESPONSE 
FUNCTIONS 

For convenience, we list here the explicit relations be- 
tween the response functions Il^and X\r%> x\> following 
from Equations (p3|). ( p4[ ) and p6| ) (omitting the super- 
script "KS"): 



n 00 
n i 
n 02 
n 03 

n X o 
n u 
n 12 
^3 

n 20 
n 2i 
n 22 
n 23 



XTT.TT + xn.u + Xll.TT + Xll, 11 
XTT.TT + XTT.1T + Xll.Tl + Xll.1T 

-*(xtt,ti _ xtt.it + xii.ti ~ xix.it) 

XTT.TT ~ XTT.U + Xll.TT ~ XlLll 
XT1.TT + XTl.ll + XlT.TT + XlT.U 

xn,n + xt Lit + xit.ti + xit.it 

-*(XT1,T1 - XT1.1T + XlT.Tl _ X1T.1t) 
XTl.TT _ XTl.ll + XiT.TT _ XlT.U 

: *(XT1,TT +XT1.U _ X1T,TT _ XlT,ll) 

: *(XTi,Ti + XT1.1T ~ XlT.Tl ~ xit.it) 

: XT1.T1 _ XT1.1T - XlT.Tl +X1T.1T 

: *(XTi,TT - XTl.U ~ XlT.TT + XlT.ll) 



1130 = XTT.TT +XTT, 11 _ Xll.TT _ Xll, U 

1131 = XTT.TI +XTT, IT _ Xll.Tl _ Xll, IT 

n 32 = -«(xtt,ti ~xtt.it _ xu,ti + xu.it) 

II33 = XTT.TT ~ XTT.il ~ Xll, TT + Xll, 11 • (Bl) 



With eigenfunctions of the form (|3l|),(p2[), the Kohn- 
Sham response function (E2h can be written as 



T',AA'( k ||^^'.^) = 5Z F ^',AA'( k lh W )^(^)^'(- 



•x.<Pi{z')ipj{z') 



(B2) 



and likewise 



n fef( k ||>2,2':W) =^(T^\\,w)ipi{z)ip j (z)ifi{z')ip j {z'). 
ij 

(B3) 

The (ki 1 , ui) are related to the F l J a , xx , (k| | , u) accord- 
ing to (Bl). The latter functions are given by 



^,AA'(k||.") 



1 ±1 



'911 



f{E s iq u ) 



(27r) 2 LO - E sl q +E s >jq I 



11] 



R*(i\\) 



\R(i\\)\ 



Vt + Oct' 4. 



i - 1 



\R(<i\\ - k n)l 

d 2 q\\ f(Esi<i U ) 



\R(d\\)\ 



l^(qn- k n)l 



, ^*(qii+ k n) 



l^(qi+ k ii)l 



V'T + Ocr'l 



\R(m)\ 



<5at + S M s 



, ^(qi+ k n) 
l^(qi+ k ii)l 



^A'T + #A'l S 



R*(m) 
\R(d\\)\ 



(B4) 



It is not difficult to show that for the case of a — > 0, i.e. vanishing Rashba term, the function F % J a , AA , becomes 



diagonal in the spin indices and reduces to the well-known limitEj (u 



Ki' AA'( k ll' w ) = -<^A<W 



d 2 g || 

(2^ 



fe?, 



i qir k n 1 ,*« , , , , ij_ qir k n 

•J 2m* m* ' L J 2m* m* 



(B5) 
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